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Abstract 

We investigate the different perturbative QCD-based models for nuclear shadow- 
ing of gluons. We show that in the kinematic region appropriate to RHIC experiment, 
all models give similar estimates for the magnitude of gluon shadowing. At scales 
relevant to LHC, there is a sizable difference between predictions of the different 
models. 



1 Introduction 



Understanding the initial stages of a ultra-relativistic heavy ion colhsion is of outmost 
importance in order to understand the outcome of the proposed heavy ion experiments, 
soon to go onhne at RHIC and later to follow at LHC Understanding the modifications 
of the parton distributions in nuclei as compared to free nucleons (shadowing) will be an 
important step towards pinning down the initial conditions of a heavy ion collision. 

At high energies (small x), there are much more gluons than any other parton species 
in a hadron/nucleus wavefunction. There are a number of processes which are sensitive to 
gluon shadowing. High energy production of mini-jets is one such example. Mini-jets will be 
important at RHIC and will dominate at LHC over soft phenomena. Nuclear shadowing of 
initial gluon distribution could significantly reduce the initial mini-jet and total transverse 
energy production. As a result, the subsequent parton thermalization will also be affected 
due to the reduced initial energy density. Production of heavy quarks is another example 
where gluon shadowing may make a dramatic difference since the probability for making a 
heavy quark pair is proportional to the square of gluon distribution function and therefore 
any depletion in number of gluons will make a significant difference in the number of heavy 
quark pairs produced. 

In recent years, there has been considerable progress made towards understanding 
gluon shadowing in perturbative QCD. Shadowing of gluons defined as 



can be understood at high energies as a recombination effect due to high gluon number 
density in the frame where the nucleus is fast, the so called Infinite Momentum Frame 
(IMF) or as a multiple scattering effect in the rest frame of the nucleus where there is a 
destructive interference between multiple scattering amplitudes. Off course, so long as one 
calculates the same physical quantity, with the same approximations made, one must get 
the same result. 

In this note, we continue our numerical study of the shadowing of gluons using 
two different QCD based formalisms; one is based on an all twist, Wilson renormalization 
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group and effective action approach to high gluon density region of QCD as developed in 
0, ^, ^. This approach takes high gluon densities into effect by including and resumming 
all n ^ 1 "hard pomeron" fusion terms in the evolution of both nuclear and nucleon gluon 
distribution functions. The other formalism is based on a generalization of the Mueller- 
Glauber multiple scattering formalism valid in the rest frame of the nucleus ^ 
(see also The two approaches lead to similar but different expressions for the gluon 
distribution function. The difference between the two approaches is investigated in |llO[ 
where it is shown that the effective action and renormalization group approach developed 
in 1^, ^, ^ is more general than the generalized Mueller- Glauber formalism f^, §|, 0; 01 
and includes effects which are not present in the latter approach. However, the difference 
between the two approaches becomes important only at very high energy (small x). Here 
we investigate the predictions of the two approaches for gluon shadowing numerically and 
show that the difference between the two is negligible in the RHIC kinematic region and 
becomes more appreciable as one goes to higher energies. For a review of the experimental 
status of nuclear shadowing, we refer the reader to . For an interesting discussion of the 
role of coherence in nuclear shadowing and its manifestation in different frames we refer 



the reader to ||12|| . 

This work is organized as follows; in the next section we briefly review the two for- 
malisms followed by a brief recall of our previous results as reported in [^. In section 2, 
we give the expressions for the nuclear gluon distribution function from the two different 
formalisms and solve the equations numerically and show our results. We finish with a 
discussion. 



2 Shadowing in IMF vs. the rest frame 

In the infinite momentum frame, shadowing can be understood as a result of high gluon 
density at a given impact parameter. The nucleus is highly Lorentz contracted due to its 
large speed. The small x gluons have large wavelengths compared to the longitudinally 
Lorentz contracted nucleus size 

1 2R 2mR 
A > — = — — 
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Therefore, small x gluons from different nucleons can spatially overlap and recombine into 
a higher x gluon. This leads to depletion of the nuclear gluon density as compared to naive 
expectation that xG^ = AxG^ . In the context of nucleons, this is referred to as saturation 
of gluon density and slows down the unlimited growth of the gluon distribution function 
which would otherwise lead to violation of the unitarity bound on physical cross sections 

In we derived an evolution equation for the gluon distribution function which 
takes — > 1 gluon ladder fusion into account and is the generalization of GLR-MQ model 



which includes only 2^1 ladder recombination. The impact parameter dependent 



evolution equation is 



r^2 2 

xG{x, Q, b±) = ^ Q'^ 



dydi 



1 A. A. 

l--exp(-)Ei(-) 

hi/ hh hh 



(2) 



where k is 



xG{x, Q, b±) 



n 3g2 

with y = log 1/x and ^ = logQ^. The exponential integral function Ei{x) is |15 

/ dt — 

Jo l + t 



(3) 



Ei(x) 



a; > 0. 



(4) 



In [m, we numerically solved this equation and calculated the nucleon and nuclear 
gluon distribution function at zero impact parameter. As written, eq. (0) is a generic 
evolution equation for gluons in either nucleons or nuclei. The distinction between nu- 
cleon and nucleus gluon distribution function is made at the initial point Xq and Qq after 
which the nucleon and nucleus gluon distribution function are determined by the evolution 
equation (see the remarks after eq. [l^). We showed in fl]] that the non-linearities in the 
evolution equation induced by the recombination effects are important. 

In the rest frame of the nucleus, shadowing is manifested through destructive interfer- 
ence between multiple scattering amplitudes as described by Glauber-Gribov type models. 
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This approach was used by Mueller to derive the following expression for the nuclear gluon 
distribution function in perturbative QCD f^, |^ 



where 



e 2 



(5) 



a^^~a;G^«^^^(x,Q2) (6) 

is the cross section for scattering of a gluon pair from a nucleon inside the nucleus. It is 
important to realize that relation holds only at the low gluon density region and will 
break down once higher twist effects become important. 

This corresponds to the following physical picture in the rest frame of the nucleus: a 
highly virtual DIS probe (or a photon) fluctuates into a gluon pair (or quark anti-quark 
pair) well before it reaches the nucleus. At small x, this pair has a long life time t ~ ^ 
and coherently scatters off the nucleons as it goes through the nucleus. The destructive 
interference between the scattering amplitudes reduces the flux of photons as seen by the 
nucleons sitting inside the nucleus which reduces the nuclear cross sections. The Mueller 
formula (^) takes into account only the interaction of the fastest (or most energetic) gluon 
pair with the nucleus. This equation was numerically solved in p, |^ and we refer the 
reader there for a comparison. It was shown in that Mueller formula leads to gluon 
shadowing which is almost independent of the initial non-perturbative shadowing input. 
It is important to realize that Mueller formula is not a non-linear equation for the gluon 
distribution function the same way GLR-MQ is. For instance, in Mueller formula, eq. 
(^) reflects a linear relation between the cross section and the gluon distribution function. 
This would not hold were there non-linear effects like higher twist terms present. It is 
also known |^ that Mueller formula over estimates the amount of shadowing. Ayala et al. 
proposed to include the effects of scatterings of the next-fastest, etc. gluon pairs with 
the nucleus by iterating the Mueller formula. To do this, one replaces the nucleon gluon 
distribution function xG^ in the exponent of eq. (§) by the nuclear gluon distribution 
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function xG^. Furthermore, assuming a Gaussian form for the shape function S{bt), they 
integrated over the impact parameter to get 0: 
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^^^^ Q) = ^ RaQ [C + H^^agl) + EiiKagl)]^ (7) 

where C ~ 0.57 is the Euler constant and 

In @, this replacement is justified by noticing that the new equation has the correct low 
density properties such that it reproduces DLA DGLAP and GLR. It should be emphasized 
that this modification of Mueller formula (|^) is just an ansatz which seems to produce cor- 
rect low density limits for the gluon distribution function and, unlike the Mueller formula, 
was not derived from QCD. It is shown in that, with a specific definition of the gluon 
distribution function, one can get equation from the F2 structure function calculated 
from the generalized Mueller-Glauber formalism. It is interesting that both approaches 
predict a slow down in the growth of the gluon distribution function at very small x such 
that 

xG{x, Q^) ~ irR^Q'^ In 1/x. (9) 

In the following, in order to compare the two equations, we will use the same Gaussian 
ansatz as used by to integrate over the impact parameter. We will then solve the two 
equations numerically starting with the same exact initial conditions and compare the 
results. 

3 Nuclear gluon distribution function 

Here we will use the Gaussian ansatz for the shape function 

Sih) = 4^e-i (10) 

SO that / (PhtS{ht) = 1 in order to perform the impact parameter integration in eq. (@). 
To do so, we first use the relation (|^) to rewrite the right hand side of eq. @ in an 
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integral form. We then switch the order of integration and perform the impact parameter 
integration first and then do the t integration. The result is 

xG{x,Q) = ^ 7rR^Q^exp{-)Ei{-) (11) 



where R is 

We now solve equations (|ll]) and (J^) numerically using the methods described in [jl]] 
in full detail. Here we will just briefly highlight the approximations made in [|l|. We 
assumed that at some initial point xq and Qo, there is no shadowing. This is motivated 
by experimental results |TT, T^] which show that the ratio ^^qn = 1 at x ~ 0.5 — 0.6 
with weak dependence. We then used the semi-classical approximation to convert our 
partial differential equations into coupled but ordinary differential equations which can 
then be solved using Runge-Kutta methods. We refer the reader to for details. The 
only difference in our choice of parameters in this work is the nuclear radius R which was 
taken to be 5/m for A = 200 in []I| while here we use a more realistic value of 7fm (more 
precisely, = RqA^^^ with Rq = 1.1 fm being the nucleon radius). To see the effect of 
integration over the impact parameter on shadowing, we first show our result for shadowing 
at zero impact parameter in Figure 1. 

The difference between the two equations as a function of x and at fixed Q = 5GeV 
is shown in Figure 2. Comparing Figures 1 and 2, it is clear that averaging over the 
impact parameter reduces the amount of shadowing as expected. Here, S'^^^^ refers to 



the shadowing ratio as defined in ([T|) as calculated from the solution of eq. ([TTD while 
gAGL calculated from the solution to eq. (0). Also, S'^^^ is shadowing calculated from 



the solution of GLR-MQ eq. [|T^ which is the second term in the expansion of eqs. ([TT]) 
and (|^. The two expressions give very similar results for RHIC {x ~ 0.01) and deviate 
appreciably only at very small x appropriate to LHC (x ~ 0.0001). For reference, we also 
show the shadowing ratio calculated from GLR-MQ. As expected, GLR-MQ predicts more 
shadowing than the other two expressions. 
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Figure 1: Shadowing from eq. 
nuclei. 



as a function of x at 64 = and Q = 5GeV for different 
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Figure 2: Model dependence of shadowing as a function of a: at Q = 5GeV. 
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We show the dependence of shadowing at fixed x as predicted by the two equations 
in Fig. 3. It is interesting to see that both approaches predict very weak dependence 
of gluon shadowing at RHIC. The difference between the two approaches becomes more 
pronounced at LHC as x gets smaller. The unphysical trend at low Q is due to both DLA 
and our semi-classical approximation breaking down at small Q as discussed in full detail 
in m. 
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Figure 3: Q dependence of shadowing ratio at RHIC [x ~ 0.01) and LHC [x ~ 0.0001). 

The dependence of our results on the choice of the parameterization chosen for the 
gluon distribution function is shown in Figure 4. In [ll], We used the CTEQ pa- 
rameterization to determine the initial gluon distribution function at scale Xq = 0.5 and 



Qo ~ IGeV. Here we use both CTEQ and GRV94 |T8| in order to compare the sensitivity 
of our results to the choice of parameterization of gluon distribution function available. 
We chose GRV94 since it is the closest in spirit to DLA approximation employed in all of 



the perturbative QCD models. Note that GRV94 is known to fail ||19[ at small values of x 
and Q but that is not relevant here since we use it at a high value of x where it is known 
to work. As seen, the two parameterizations predict similar results for RHIC kinematic 
region while for LHC there is a bigger difference. 
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Figure 4: parameterization dependence of shadowing ratio from eq. (|TT1) . 

4 Discussion 

We discussed the different perturbative QCD-based models for nuclear gluon distribution 
function and numerically investigated their predictions for x and Q dependence of shad- 
owing in the kinematic region appropriate to RHIC and LHC We showed that predictions 
of different models for shadowing of gluons at RHIC are comparable while the difference 
at LHC can be of order ~ 10% for Gold or Lead. An important point which needs to 
be investigated further is inclusion of initial shadowing of gluons due to non-perturbative 
effects at the starting point of evolution xq after which the perturbative evolution takes 
place. This is currently under investigation and will be reported on later PU|. The more 
experimentally relevant quantity to investigate is the shadowing of the nuclear structure 
function F.f since shadowing of gluons is not directly observable. The all twist and 
as well as the longitudinal structure function Fl were computed in at the classical 
level. Including the quantum loop effects due to gluons is straight forward and is currently 



under investigation |2T]]. One can then predict the experimentally measured shadowing 
ratio F^/F^ for different nuclei at different x, as well as the longitudinal structure 
functions. 
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